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Abstract 

By using the most sensitive two-point correlation functions introduced to date, we reconstruct 
the microstructures of two-phase random media with heretofore unattained accuracy. Such media 
arise in a host of contexts, including porous and composite media, ecological structures, biological 
media, and astrophysical structures. The aforementioned correlation functions are special cases 
of the so-called canonical n-point correlation function H n and generalize the ones that have been 



q | recently employed to advance our ability to reconstruct complex microstructures [Y. Jiao, F. H. 



Stillinger, and S. Torquato, Proc. Nat. Acad. Sci. 106, 17634 (2009)]. The use of these gen- 
eralized correlation functions is tantamount to dilating or eroding a reference phase of the target 
medium and incorporating the additional topological information of the modified media, thereby 
providing more accurate reconstructions of percolating, filamentary, and other topologically com- 
plex microstructures. We apply our methods to a multiply-connected "donut" medium and a 
dilute distribution of "cracks" (a set of essentially zero measure) demonstrating improved accuracy 
in both cases with implications for higher-dimensional and biconnected two-phase systems. The 



X 

■ high information content of the generalized two-point correlation functions suggests that it would 

be profitable to explore their use to characterize the structural and physical properties of not only 
random media, but also molecular systems, including structural glasses. 
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I. INTRODUCTION 



Two-phase random media, which are partitions of space into two distinct phases with 
interfaces that are known only probabilistically, abound in a host of situations, including 



composite materials jl-5], cosmological (6|, |7| and geophysical structures and biological 
systems jg, 9]. Importantly, the effective physical and chemical properties of random media 
are rigorously related to microstructural correlations functions |2j. In principle, complete 
microstructural information requires an enormously large set of such correlation functions, 
including statistical descriptors describing higher-order correlations (e.g., three-point, four- 
point, and higher-order statistics). Since in any practical setting only lower-order correlation 
functions are available, it is of paramount interest to quantify the amount of microstructural 
information contained in such reduced statistics and then to identify the most sensitive 
descriptors. 

Microstructure reconstruction techniques provide a means of probing these issues. Re- 
construction of a random medium that matches limited structural information, i.e., finite 



set o 
lem 



targeted lower-order correlation functions of the system, is an intriguing inverse prob- 
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11] . By comparing unconstrained correlation functions between the original and 



reconstructed medium, one can quantify the capacity of the reduced statistics to specify 
the important structural features of the medium 12|. An effective reconstruction proce- 
dure enables one to generate accurate renditions of random media at will, and subsequent 
analysis can be performed on the reconstruction to obtain desired macroscopic properties 
nondestructively 

Although there has been significant progress on reconstruction algorithms 0, Q-Q, 



it remains a challenge to reconstruct accurately percolating microstructures, sets of essen- 
tially zero measure (e.g., cracked media and filamentary structures), and multiply connected 
microstructures. In this paper, we introduce a procedure that enables one to begin to re- 
construct such microstructures with heretofore unattained accuracy. 

It is noteworthy that a panoply of different types of correlation functions arise in rig- 
orous theories of structure-property relations for heterogeneous media, which includes the 
n-point probability function S n , n-point surface correlation functions, pore-size functions, 
lineal-path functions, nearest-neighbor functions, among others [l|. It was shown that all 
of these statistical descriptors are special cases of the more general canonical n-point cor- 



relation function H n , which is a hybrid correlation-density function and can be represented 
analytically for a wide class of microstructures . Indeed, H n characterizes a modified 

random medium that involves uniformly "dilating" or "eroding" a reference phase in the 
normal direction to the two-phase interface in the actual random medium Jj. Dilation is 
achieved by placing a "test" sphere of radius 5 into the non-reference phase and determin- 
ing the space available to the sphere as if the reference phase was impenetrable to the test 
sphere. The space unavailable to the test sphere is considered to be the modified or dilated 
reference phase. The erosion process simply reverses this procedure such that the reference 
phase is allowed to be penetrated by a distance 5 normal to the two-phase interface. The 
quantity H n is an even more general quantity since it allows for multiple test spheres, each of 
which has a different size; see Refs. |l| and 17| for details. Thus, H n contains considerably 



more information than any of the types of correlations functions that arise in the theory 
of heterogeneous media. We will use special cases of H 2 in our reconstruction method, as 
described below. 

In practice, the enormous information content of the full n-point canonical function H n 
is neither experimentally nor theoretically accessible, and one must instead rely on reduced 
statistical information about the microstructure. It is an open and fundamental problem 
to identify via reconstruction techniques lower-order statistical descriptors that can be both 
manageably measured and yet reflect a substantial portion of the complete microstructural 
information contained in H n . The most prominent example of reduced statistical information 
in the two-point probability function ^(r), probability of finding two points separated by 
the displacement vector r in one of the phases. This function is accessible from scattering 



experiments 



18] or tomography [19[. Such information alone has been shown to be generally 
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12| 



insufficient to render accurate reconstructions 

The most successful reconstruction procedures to date have directly incorporated con- 
nectedness information about the target media. In particular, Jiao, Stillinger and Torquato 



15[ have shown that incorporation of the two-point cluster function C 2 , experimentally 
obtainable from tomography or other three-dimensional imaging techniques results in 
superior reconstructions. The two-point cluster function defines the probability of finding 
-wo points in both the same phase and in the same connected cluster of the microstructure 



1. 



l|, |20[. The function C 2 {v) is the connectedness contribution to the two-point probability 
function, namely, S 2 {r) = C 2 (r) + E 2 (r), where E 2 is the two-point blocking function, cor- 



responding to the case where two points fall within different clusters of the same phase. It 
has been shown that C2 contains substantially more information in excess to S2 than even 
the three-point probability function S3, thereby highlighting the inefficiency of incorporating 



15|- 



additional information via the usual higher-order correlation functions 

However, there remains a large class of microstructures that are difficult to characterize 
accurately using only S 2 and C 2 . Examples of such systems include heterogeneous media 
with percolating and/or multi-connected phases 21], systems with zero-measure structures 
such as cracked materials and filamentary structures that arise in the large-scale structure of 



the Universe 



221 ] . and other topologically complex media. In particular, C2 works best when 



the phase of interest is below its percolation threshold; otherwise, it is a long-ranged function 
that does not reflect structural details. Moreover, when the phases occupy effectively a space 
of Lebesgue measure zero, the probability of finding a point "within" the set identically 
vanishes; these elements are therefore unable to contribute to either S2 or Ci- 

In this paper, we present an improved reconstruction procedure that naturally generalizes 
the S2.-C2 methodology presented in Ref. [l||. Fundamental to our work is the informa- 
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241 ] . In particular, 



tion contained in the void space external to the phase of interest 
we incorporate special and new cases of H 2 , which contain additional topological informa- 
tion associated with modifications of the target microstructure via "dilation" and "erosion" 
processes, thereby allowing a more accurate reconstruction procedure. Importantly, our 
methods contain previously-reported results for S2 and C2 as a special case, thereby provid- 
ing a general means to study the amount of information contained in lower-order correlation 
functions. 



II. GENERALIZED TWO-POINT CORRELATION FUNCTIONS: SPECIAL CASES 
OF H 2 



Previous studies on reconstructions of heterogeneous media 15| have utilized the pore-size 
probability function F(S), defined as the probability of finding a spherical cavity of radius 
5 contained in the void space external to the phase of interest. The pore-size probability 
function contains short-range connectedness information within the void space and can also 
be interpreted as the volume fraction of the pore space after uniformly dilating the reference 



phase by a linear scale 5 



23]; note that erosion by a linear scale 5 corresponds to the 



continuous extension of F to negative values of 8. As such, it is a sensitive one-point 
descriptor of heterogeneous media implicitly containing limited higher-order information. 

Higher-order pore-size probability functions are simply special cases of the canonical 
ra-point correlation function H n . In particular, the n-point pore-size probability function 
F n (r n ; {8} n ) is the probability of finding n spherical cavities of radii 8 
the "void" space external to the reference phase; in the notation of Ref. 
if„(r n ; 0; 0). Of particular interest here is the two-point pore-size probability function 
F 2 (r; 8) for a statistically homogeneous medium, where r = r 2 — r x and 8 — 61 — 82- This 
function also contains short-range connectedness information about the modified void space 
and can be interpreted as the two-point probability function S2 for the void space after 
uniformly dilating the reference phase by a linear scale 8. As such, it can be decomposed 
as F 2 (r; 8) = K 2 (r; 8) + M 2 (r; 8), where K 2 and M 2 are the two-point pore-size cluster and 
blocking functions, respectively. 

Since F 2 (r;8) — > 5*2 (r) (for the void space) as 8 — > 0, it follows that F 2 and its cluster 
counterpart K 2 rigorously contain more information about the microstructure than the stan- 
dard £2 and C 2 functions alone. Thus, we incorporate these generalized S2 and C 2 functions 
in reconstructions, which amounts to applying dilation or erosion processes to modify the 
phase connectedness of a target microstructure. We then perform a reconstruction proce- 
dure on the modified medium, which is expected to yield improved reconstructions of the 
original unmodified medium. 
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III. RECONSTRUCTION PROCEDURE AND APPLICATIONS 



Our reconstruction procedure is a modification of the Jiao-Stillinger-Torquato methodol 



ogy [15|, which extends the Yeong-Torquato simulated annealing algorithm to include 
clustering information. In particular, we overlay the image with a set of N d pixels taking 
values 1 (within the phase of interest) or (in the void space), where d is the Euclidean di- 
mension. By measuring the distances between void and surface pixels, we perform an initial 
dilation or erosion process to modify the microstructure prior to reconstruction. Recon- 
struction of the modified structure then proceeds by stochastic optimization of the following 




a 



(b) 



(c) 



FIG. 1. (Color online) The doubly-connected "donut" medium, (a) Original microstructure (outer 
rings) with modified target structure (inner disks; red online) after an erosion process, (b) Direct 
S2(r)-C2(r) reconstruction of the original microstructure. (c) 5 , 2(r)-C l 2(r) reconstruction of the 
modified target structure (inner disks; blue online) along with the final reconstructed system (outer 
rings) after a dilation process. All reconstructions have a final energy (error) at most E ~ 1CF 8 . 

objective (energy) function 

M S 
a=l j=l 

where M is the number of sampling points and 5 is the number of constrained correlation 
functions fj with target values /j. In this work we consider radially-averaged functions [e.g., 
52 (r)], which are more readily accessible from scattering experiments and more relevant to 
studies in the thermodynamic limit 12[ . However, our methods can be directly extended to 
include anisotropic correlation functions. The Jiao-Stillinger-Torquato procedure tracks the 
occurrence of "cluster" and "surface" events to update efficiently the cluster and boundary 
lists at each step of the optimization procedure [15] . A cooling schedule is chosen such that 



25]. 



the energy E approaches its ground-state value within a very small tolerance 

We illustrate our improved reconstruction methodology by studying the doubly-connected 
"donut" structure shown in Figure [U consisting of annuli distributed in the plane. It is note- 
worthy that nanoparticles of mesoporous silica shells with hollow interiors have been designed 
with similar cross- sectional structures for applications to biomedical imaging and targeted 
drug delivery 26]. The uniformity of the nanoparticle sizes and shapes, including the hollow 
interiors, is particularly important and can be directly quantified using a combination of 



imaging techniques and pore-size distribution functions 



26] 




(a) (b) 

FIG. 2. (Color online) (a) S^r) and C2(r) for the modified target disk structure following an 
erosion process and for the reconstructed disk microstructure. (b) Lineal path functions L(z) for 
the original target donut structure, a direct 5 , 2(r)-C , 2(r) reconstruction of the donut system, and 
the reconstructed material using erosion and dilation processes. 

The connectedness information of the donut structures is nontrivial because the two- 
dimensional annulus is not a simply-connected topological object. A direct reconstruction 
using S2 and C2 shows that, while the two-point cluster function is able to identify sit- 
uations where points fall within the same cluster, it cannot distinguish situations where 
the "donut shells" are closed and open. Furthermore, isolated disks of particles appear in 
the reconstructed image with diameter approximately equal to the length over which C2 is 
nonzero. 

We are able to improve the reconstruction procedure by targeting the spatial information 
contained in the "cores" of the ring structures. Although the annuli surrounding the cores 
are not simply-connected, by performing an erosion process on the initial microstructure, 
we generate a "modified" system composed of simply-connected, spatially separated disks in 
the plane (see Figured]). Reconstruction of the modified medium using S2 and C2 is able to 
reproduce the statistical distribution of the particle centers, and after applying a subsequent 
dilation process, the original microstructure is recovered with greater fidelity (Figures [1] and 

ED. 

We quantify the accuracy of the reconstruction procedure by measuring some uncon- 
strained correlation function not targeted by the reconstruction algorithm. In Figure [2j we 
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compare the lineal path function L(z), defined as the probability of finding a connected 
linear segment of length z contained completely in the reference phase, for the original 
donut structure, the direct S2-C2 reconstruction, and our reconstruction using erosion and 
dilation processes. The direct reconstruction algorithm is not able to capture accurately 
this type of short-range continuous connectedness information, accounting for the poorer 
reconstruction in Figure [TJ However, the lineal path functions of the original system and 
the reconstructed material utilizing erosion and dilation processes match closely, implying 
that this information is already implicitly encoded in the spatial distribution of the donut 
centers. By utilizing effective two-point correlation functions that directly probe the cores 
of the donut inclusions, we are able to reduce the topological complexity associated with the 
reconstruction procedure without losing fundamental information about the connectedness 
properties of the microstructure. 
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(a) (b) (c) 

FIG. 3. (Color online) (a) Target microstructure consisting of a dilute distribution of cracks of 
nearly-zero measure, (b) Direct Si{r)-Ci{r) reconstruction of the target cracked structure, (c) 
Si{r)-Ci(r) reconstruction of the dilated target microstructure followed by an erosion process. All 
reconstructions have a final energy (error) at most E ~ 10 -8 . 

We have also applied our methodology to the reconstruction of a dilute distribution of 
"cracks" of nearly- zero measure, shown in Figure [31 Note that this prototypical cracked 
material is very difficult to reconstruct directly, even utilizing the two-point cluster function 
(see Figure [3b)- Although the cracks are given a finite size by the pixelation procedure, 
the volume fraction associated with the cracks is small, and a direct S2-C2 reconstruction 
preserves neither the cracked elements of the target structure nor the void-space distribution. 

However, by dilating the target structure and reconstructing the dilated system, we have 
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(a) (b) 

FIG. 4. (Color online) (a) ^(r) and C^?") for the dilated target cracked structure and its recon- 
struction, (b) Lineal path functions for the dilated target cracked structure, its reconstruction, 
and the dilation of the direct S2{r)-C2{r) reconstruction of the initial target microstructure. 

been able to improve the accuracy of the reconstruction procedure. Figure H] shows excellent 
agreement between the two-point probability and cluster functions for the dilated target and 
reconstructed microstructures, and calculation of the lineal path functions demonstrates that 
the dilation methodology provides quantitatively substantial improvement in the accuracy of 
the reconstruction. As previously mentioned, extended filamentary structures arise in studies 
of the large-scale structure of the Universe p,0|; an< ^ our results therefore suggest that two- 
point pore-size correlation functions contain important structural information about such 
systems, beyond even the standard two-point probability function S 2 and the two-point 
cluster function C 2 . 



IV. DISCUSSION AND CONCLUSIONS 



We have introduced a different reconstruction procedure that utilizes two-point pore-size 
correlation functions embodied in if 2 in conjunction with clustering information to improve 
reconstructions of two-phase random heterogeneous media. It is noteworthy that these two- 
point functions and the canonical n-point function H n are generalizations of the Minkowski 
functionals 
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. Our results support the important notion that the voic 



to the microstructural properties of a heterogeneous system 
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space is fundamental 



24J. Our methodology 



provides a natural extension of the Jiao-Stillinger-Torquato reconstruction algorithm and 



the correlation functions used therein 15]. A major conclusion of this earlier study was 
that the two-point cluster function C2 employed actually embodies higher-order structural 
information in a way that makes it a highly sensitive statistical descriptor over and above 
the standard two-point function S^. Thus, C2 has the ability to leapfrog past the usual 
approach of incorporating additional information via higher-order versions of S2 (i.e, S3, 



5*4, etc.) 29]. The present work clearly distinguishes the two-point probability and cluster 
functions F2 and K2, respectively, as highly sensitive descriptors of random media, since 
they are generalizations of S2 and C2 and hence contain even more information than the 
combination of S 2 and C 2 . 

The high information content of F 2 and K 2 suggests that it would be profitable to explore 
their use to characterize the structure and physical properties of not only random media, 
but molecular systems. For example, these functions might provide sensitive structural 



signatures of glassy states of matter (an issue of current interest 30]) and improved rigorous 
upper and lower bounds on the effective properties of a heterogeneous material jl|. These 
and related issues are beyond the scope of the present work, and hence represent interesting 
future research topics. 

As noted earlier, two-phase media with percolating phases are difficult to target directly 
with C2, which becomes long-ranged at the percolation threshold. However, preliminary 
reconstruction results for such cases obtained by applying dilation and erosion processes to 
generate a nonpercolating target medium, which is then reconstructed using the generalized 
5 2 and C 2 , are promising and will be discussed further in future work. 

Although we have focused in this paper on reconstructions of two-dimensional heteroge- 
neous media in order to justify our methods, previous work on reconstructions using S2 and 



C2 by Jiao, Stillinger, and Torquato [15] suggests that our methods can be easily extended 
to higher dimensions. One interesting three-dimensional application of our work is the re- 
construction of biconnected media, in which both phases percolate simultaneously. This 
behavior, though difficult to realize in two dimensions, is more common in three dimensions 
1,|10J, with important implications for the effective properties of a heterogeneous material 
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